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ABSTRACT 

Hot Jupiter atmospheres exhibit fast, weakly-ionized winds. The interaction of these winds with the planetary 
magnetic field generates drag on the winds and leads to ohmic dissipation of the induced electric currents. We 
study the magnitude of ohmic dissipation in representative, three-dimensional atmospheric circulation models 
of the hot Jupiter HD 209458b. We find that ohmic dissipation can reach or exceed 1% of the stellar insolation 
power in the deepest atmospheric layers, in models with and without dragged winds. Such power, dissipated in 
the deep atmosphere, appears sufficient to slow down planetary contraction and explain the typically inflated 
radii of hot Jupiters. This atmospheric scenario does not require a top insulating layer or radial currents that 
penetrate deep in the planetary interior. Circulation in the deepest atmospheric layers may actually be driven 
by spatially non-uniform ohmic dissipation. A consistent treatment of magnetic drag and ohmic dissipation 
is required to further elucidate the consequences of magnetic effects for the atmospheres and the contracting 
interiors of hot Jupiters. 
Subject headings: 



1. INTRODUCTION 

The discovery of the first transiting extrasolar planet 
HD 209458b (Charbonneau et al. 2000) has opened a new 
chapter for the study of planetary bodies. This planet belongs 
to the class of hot Jupiters, which are close-in gaseous giant 
planets thought to be tidally locked to their host star. In the 
last decade, a wealth of observations has allowed direct in- 
vestigations of the atmospheric and bulk properties of these 
planets, while much theoretical work has been aimed at the 
interpretation of these data (see, e.g., Showman et al. 2010, 
Burrows & Orton 2010, Baraffe et al. 2010 and Seager & 
Deming 2010 for recent reviews). A property of hot Jupiters 
which has constituted a longstanding puzzle is their anoma- 
lously large radii, which have been interpreted as requiring 
that extra heat be deposited in the convective regions, or alter- 
natively in the deep atmospheres (at pressure of tens of bars, 
e.g. Guillot & Showman 2002) of these planets. 

In recent years, progress in modeling the unusual atmo- 
spheric circulation regime of hot Jupiters, with permanent day 
and night sides, has also been made (e.g. Cooper & Show- 
man 2005; Dobbs-Dixon & Lin 2008; Showman et al. 2009; 
Rauscher & Menou 2010; see Showman et al. 2010 for a 
review). These purely hydrodynamical models have gener- 
ally ignored the possibility that magnetic effects acting on the 
weakly-ionized winds could significantly influence the circu- 
lation pattern. 

Recently, Perna et al. (2010) evaluated the level of mag- 
netic drag on a representative hot Jupiter atmospheric flow, 
and argued that it is likely to provide an effective frictional 
mechanism to limit the asymptotic speed of winds in these 
atmospheres. Batygin & Stevenson (2010) studied the role 
of ohmic dissipation in hot Jupiters, using a prescribed zonal 
wind profile and a one-dimensional atmospheric structure 
model, to argue that the typical amount of heat deposited deep 
in these planets' convective zones is sufficient to explain their 
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inflated radii. In this Letter, we build on our previous work 
on magnetic drag (Perna et al. 2010; hereafter PMR10) to 
investigate the magnitude of ohmic dissipation in the atmo- 
spheres of hot Jupiters, and its consequences for the dynam- 
ics and the thermal evolution of these planets. We use spe- 
cific three-dimensional atmospheric circulation models of the 
planet HD 209458b. We compare the amount of Ohmic heat- 
ing expected for typical magnetic field strengths with the extra 
heat required in the deep atmosphere to slow down contrac- 
tion according to planetary evolutionary models (Guillot & 
Showman 2002). While our calculations are specific to the 
case of HD 209458b, our overall results are expected to hold 
more generally for hot Jupiters with similar gravity, irradia- 
tion strength, and magnetic field. 

2. ATMOSPHERIC CURRENTS AND OHMIC DISSIPATION 
2.1. Atmospheric models and ionization balance 

The fiducial atmospheric circulation model used here was 
computed by Rauscher & Menou (2010) for HD 209458b, 
under the assumption of no significant magnetic drag on the 
atmospheric flow. The model describes the atmospheric flow 
in a frame that is rotating with the bulk planetary interior. The 
meridional and zonal wind speeds in the atmosphere, as well 
as its thermodynamic variables, are returned at each grid lo- 
cation in the three-dimensional model atmosphere. Location 
is identified by the angular spherical coordinates (9,<j>) and 
pressure, p, for the vertical coordinate, The model bottom 
is located at 220 bar while the top level is set at a pressure 
of 1 mbar. In addition to this drag-free model, we also per- 
form some of our calculations for the model with strongest 
drag described in PMR10. Since, apart from wind drag, these 
two models are identical, this allows us to evaluate the con- 
sequences for ohmic dissipation of an atmospheric flow with 
significantly dragged winds. 

In our circulation models, the local heating/cooling rate (en- 
ergy per unit mass) is modeled as Newtonian (linear) relax- 
ation, Q T = (T eq (p, 6», (/))-T)/T md (p), where r rad represents the 
radiative timescale on which the local temperature T relaxes 
to the prescribed equilibrium profile r eq . The nature of the 
atmospheric circulation obtained with this simplified forcing 
compares well to what is obtained with more realistic forcing 
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(e.g., Showman et al. 2009). Following Cooper & Showman 
(2005), Rauscher & Menou (2010) relied on the work of Iro 
et al. (2005) to implement the detailed profiles of this radia- 
tive forcing. The atmosphere is divided into actively forced 
layers, above the 10 bar level, and "inert" layers below. The 
radiative forcing is assumed to be negligible in the inert lay- 
ers, which corresponds to r la d — > oo. The active layers, on the 
other hand, are forced on a finite radiative timescale, r la d(p). 
The local atmospheric temperatures obtained after model re- 
laxation hover at about 1800 K in the deepest levels, while in 
the upper levels they range from about 500 K on the night side 
to about 1500 K on the day side (see Rauscher & Menou 2010 
for details 4 ). 

At these temperatures, the primary source of free electrons 
stems from thermally ionized alkali metals with low first ion- 
ization potentials: Na, Al, K. For simplicity here (and con- 
sistently with PMR10), we adopt an approximation to Saha's 
equation (Balbus & Hawley 2000) which assumes potassium 
to be the dominant contributing species: 
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Here n e and n„ are respectively the electron and neutral num- 
ber densities (in cm" 3 ) and a K ~ 10~ 7 is the potassium solar 
abundance. As discussed in PMR10, equation (JTJ is a good 
approximation as long the resulting ionization fraction, x e , is 
<C cik\ this condition is satisfied in our atmosphere models. 

We assume that the gas is overall neutral, which implies an 
equality between the electron number density n e and the ionic 
one The electrical conductivity and associated resistivity 
are given by (see also Laine et al. 2008) 
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respectively, with the collision rate between electrons and 
neutrals approximated as (Draine et al 1983) 
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2.2. Induced currents and Ohmic dissipation 

The first step towards the computation of currents involves 
an estimate of the importance of various non-ideal MHD 
terms in the full induction equation for the weakly-ionized 
medium. For the atmospheric models under consideration 
here, PMR10 found that, for a surface magnetic field strength 
of 3 G, the Hall term is completely negligible throughout 
the flow, while the resistive term largely dominates over the 
ambipolar diffusion term everywhere except possibly in the 
model uppermost levels. Therefore, to a good approximation, 
the induction equation can be considered as purely resistive. 

4 While the simulations by Rauscher & Menou (2010) adopt an approx- 
imate treatment for radiative transfer, other simulations by Showman et al. 
(2009) with explicit radiative transfer find rather similar results for the day- 
night temperature gradient and wind pattern under comparable physical con- 
ditions. However, it should be recognized that all these models are subject to 
some uncertainties for the development of the deep atmospheric winds due to 
the very long integration times needed for spin-up. 
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FIG. 1 . — Cylindrical maps of Joule heating times (logry, in sec), at four 
pressure levels in our fiducial, drag-free atmospheric model. The sub-stellar 
point is centered at longitude and latitude zero. From top to bottom, Joule 
heating times are shown at 1 mbar, 50 mbar, 2 bar, 90 bar. Ohmic dissipation 
is not spatially uniform and it can dominate over heating by stellar insolation 
at deep enough pressure levels (p > 2 bars), where a new form of Joule- 
driven circulation could emerge. The calculation was performed for B = 3 G, 
but note that ohmic dissipation increases steeply with magnetic field strength 
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Assuming that zonal winds are dominant, one needs only to 
consider the toroidal component of the induction equation 
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where = v ( p{rsm9)~ l is the local angular velocity of the 
flow, r is the radial spherical coordinate, and the magnetic 
field is assumed to be an axisymmetric dipole. 

In PMR10, we also argued that, to leading order, the lati- 
tudinal component of the current induced by the zonal flow 
should be dominant. This assumption could be verified in 
more general versions of our models. Since the magnetic 
Reynolds number is R m <C 1 for the bulk of the flow, no sig- 
nificant dynamo action is expected in the atmosphere. The 
dipolar planetary field {B r ,Bg) is maintained by currents in 
the interior of the planet, while the zonal flows separately in- 
duce a component from the dipolar one in the superficial 
atmospheric layers. Under these conditions (Liu et al. 2008), 
the resulting steady-state latitudinal current can be computed 
as: 

, „ , s csin6> f R , , ,2 fan i ao. 
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where the last term includes a boundary current in the upper- 
most modeled level, jg(R,6,<f>). Lacking information about 
the nature of currents possibly flowing from regions above the 
modeled atmospheric layers, we set this boundary current to 
zero for simplicity. As discussed in PMR10, this unknown 
boundary current represents an important source of uncer- 
tainty in our modeling, but, unless near cancellations occur, 
additional boundary currents could in principle contribute to 
even stronger ohmic dissipation than estimated here. 

The ohmic power per unit volume dissipated locally is read- 
ily computed as (e.g., Liu et al. 2008) 
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2.3. Joule-driven Circulation 



We first examine the possibility that spatially non-uniform 
ohmic dissipation drives a circulation in regions of the atmo- 
sphere where it is comparable or stronger than heating due 
to stellar insolation. To permit a direct comparison between 
the two heating sources, we find it convenient to define a typ- 
ical local timescale associated with Joule heating, which is 
deduced from the energy equation 
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so that the Joule heating time is 
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The specific heat of the gas in our model atmospheres is C p = 
1.43 x 10 8 erg g" 1 K" 1 . 

Fig. 1 shows cylindrical maps of Joule heating times for 
a nominal surface dipolar field of 3 G, at four depths in the 



atmospheric flow of our fiducial, drag-free model. In the up- 
permost levels, the Joule heating times span almost 20 orders 
of magnitude, with the shortest times found on the day side. 
This wide span reflects the large variations in resistivity, r], 
between the day and the night sides. On the day side, which 
has higher temperatures, rj is considerably smaller than on 
the night side, yielding larger currents, which in turn result 
in larger Qj and correspondingly smaller tj values. On the 
other hand, deep in the inert layers, which are little affected 
by stellar irradiation, Joule heating times span a more mod- 
est range of values, ~ 10 9 — 10 12 s over a large fraction of the 
flow. Longer timescales along the equator in the deeper levels 
are mostly the result of geometric effects: in Eq.d5]l, the term 
(dfl/dr)B r generally dominates over 1 /r(dQ/d9)Bg. Since 
Q oc sin -1 9 and B r oc cos 9, the net result is j oc cos 9, and 
hence tj oc cos" 2 #. For a non-axisymmetric dipolar field, the 
anisotropy pattern would be different. 

To evaluate the possibility of Joule-driven circulation, we 
compare local Joule heating times with a representative heat- 
ing time associated with stellar irradiation, T m . For con- 
sistency with the Newtonian forcing scheme used in our 
atmospheric models, we adopt a simple downward insola- 
tion flux that approximately matches the absorption prop- 
erties of the one-dimensional models computed by Iro et 
al. (2005) at a few bars level. The flux is taken to obey 
F{p) =F exp[-4.5/V(2bar)], where F = 2.23 x 10 8 erg cm -2 
s _1 is the incident flux at the model top for HD 209458b, after 
geometric dilution by a factor 1/4 as is customary for ID- 
averaged models (e.g., Hansen 2008). A stellar irradiation 
heating timescale, r m , is then computed following Eq. ([8]), 
with a heating rate taken as the vertical divergence of the stel- 
lar irradiation flux, Q\ n = -dF(r) /dr. This allows us to eval- 
uate r m simply as a function of the local pressure, pressure 
scale height and density. For example, at p = 2 bar, using 
p ~ 5.5 x 10" 5 g cm" 3 , we estimate r ilT ~ 5 x 10 7 s. At p = 10 
bar, where the stellar insolation flux has been very strongly at- 
tenuated already (Iro et al. 2005), we estimate T m ~6x 10 15 s 
for p ~ 3 x 10~ 4 g cm -3 . Deeper in the atmosphere, the inso- 
lation flux is further attenuated (exponentially so) and heating 
by insolation becomes largely inconsequential. 

A comparison between values for r m and the detailed maps 
of tj values in Fig. 1 shows that, for the nominal 8 = 3G used 
in our tj calculation, irradiation dominates over Joule heating 
(7"irr <C Tj) everywhere on the day side of the atmosphere at 
pressure levels above a few bars. In these upper regions, we 
expect the atmospheric thermal structure to be largely unaf- 
fected by the extra Ohmic heating. Deeper than a few bars, 
however, Joule heating times start becoming comparable to 
or shorter than the typical heating time associated with stellar 
insolation, over significant regions of the atmospheric flow. 
Eventually, at levels deeper than about 10 bar, ohmic dissipa- 
tion easily dominates over stellar insolation over much of the 
atmosphere. Since this extra source of heating in the deeper 
regions of the atmosphere is spatially non-uniform (Fig. 1), 
it should lead to some form of Joule-driven circulation in 
the "inert" layers located well below the radiatively-forced 
"weather" layers. This conclusion is largely independent of 
the specific stellar insolation model adopted here since the 
existence of radiatively inactive layers at levels deeper than a 
few bars is a rather generic property of irradiated hot Jupiter 
atmospheres (e.g., Hansen 2008). However, it remains to be 
seen what the nature of such Joule-driven circulation might 
be given the uneven heating pattern shown in Fig. 1, the non- 
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uniformity of stellar insolation from the day-side to the night- 
side and the presence of an additional, spatially uniform net 
flux emerging from the deep planetary interior. Finally, it is 
worth remembering that, since tj oc B~ 2 , a stronger field could 
produce a substantially more dominant ohmic dissipation than 
estimated above for B = 3 G 5 . 

24. Inflating Hot Jupiters 

Another potentially important consequence of ohmic dis- 
sipation is the possibility that it contributes to the tendency 
for hot Jupiters to have inflated radii (Batygin & Stevenson 
2010). To evaluate the magnitude of this effect on the basis of 
our three-dimensional atmospheric models, we compute the 
cumulative ohmic power dissipated in successive atmospheric 
shells, from the top level at a pressure p, of 1 mbar, down to 
a pressure p(r) within the atmosphere; this is readily obtained 
by integration of Eq.([6]i, 
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Note that the power as computed above includes a horizontal 
average over the day and night sides. This averaging washes 
out the horizontal features in Joule heating shown in Fig. 1, an 
issue deserving further attention in future multi-dimensional 
models. Fig. 2 displays this ohmic power as a function of 
pressure for our fiducial, drag-free atmospheric model and for 
the model with strongest drag described in PMR10, with zonal 
wind speeds typically reduced by about 30-50%. For each 
model, we display ohmic power for two different values of 
the magnetic field strength adopted in the ohmic dissipation 
calculation, 8=3 G, and B = 10 G. This allows us to sepa- 
rately explore the effects of drag on the zonal winds and of the 
magnetic field strength on the resulting magnitude of ohmic 
dissipation. This is a useful exercise until models with a self- 
consistent treatment of magnetic drag and ohmic dissipation 
become available. 

Fig. 2 shows that much of the dissipated ohmic power 
builds up at pressure levels from a few bars to several tens 
of bars. The power scales as B 2 and, relative to the drag-free 
model, it is reduced by a factor ~ 3 in the circulation model 
with strongly dragged winds, at fixed magnetic field strength. 
Since the stellar insolation power received by HD 209458b is 
~ 3 x 10 22 W, Fig. 2 shows that the total ohmic power dis- 
sipated in these model atmospheres typically approaches or 
exceeds 1 % of the stellar insolation power. Notice that, while 
the Ohmic power in Fig. 2 is integrated down to the lowest 
grid points (~ 200 bar) available in the model of Rauscher 
& Menou (2010), it is not yet seen to saturate at those lev- 
els. This is because, while the zonal velocities decrease with 
depth, they are still non-zero in the deepest, 'inert' model lay- 
ers. It is in fact not a priori clear where exactly the transition 
to the convection zone should be located in hot Jupiters and 
this constitutes an uncertainty for Joule heating models of the 
deepest atmospheric layers. 

Overall, our results on the dissipated ohmic power are 
consistent with the estimates made by Batygin & Stevenson 

5 The magnetic field strengths of hot Jupiters are unconstrained from an 
observational point of view. In the case of Jupiter, the measured field is 14 G 
at the pole and 4.2 G around the equator. Arguments (e.g. Christensen et al. 
2009) suggest that the field strength should scale with the square root of the 
density and the spin frequency, which could make the B field in hot Jupiters 
somewhat weaker. Calculations of the internal structure and convective mo- 
tions of giant planets (Sanchez-Lavega 2004) yield surface magnetic fields 
~ 1 — 2 G for synchronized hot Jupiters. 
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FIG. 2. — Cumulative ohmic power dissipated above pressure level p, in 
four different models for HD 209458b. The solid lines correspond to the drag- 
free model described in Rauscher & Menou (2010) while the dashed lines 
correspond to the model with strongest drag, and reduced winds, described 
in Perna et al. (2010). Curves are labeled with the value of the magnetic field 
adopted for the ohmic dissipation calculation. The ohmic power can reach 
or exceed 1% of the stellar insolation power, ~ 3 X 10 22 W, in the deepest 
atmospheric layers, in models with and without dragged winds. 

(2010) on the basis of a more idealized, one-dimensional at- 
mosphere model, using a parametrized zonal wind profile. 
Like them, we find that a large fraction of the ohmic power is 
dissipated in the deeper atmospheric levels, below ^10 bar. 
Contrary to Batygin & Stevenson (2010), however, we do not 
solve for radial currents in our models, only meridional ones 
(Eq. and our results are thus largely independent of the 
presence or absence of an insulating layer high in the atmo- 
sphere. In fact, we find that an insulating layer is present on 
the nightside in our three-dimensional models, but not neces- 
sarily on the much hotter, less resistive dayside. This raises 
the possibility that, in more detailed three-dimensional MHD 
models, current loops could actually close high in the atmo- 
sphere, rather than deep in the planetary interior as conjec- 
tured by Batygin & Stevenson (2010). While our results, re- 
lying on leading-order meridional currents in the atmospheric 
region alone, are presumably not strongly sensitive to these 
conditions affecting radial currents, it remains to be seen how 
currents would flow in more realistic non-axisymmetric and 
three-dimensional models. 

Various authors have evaluated the extra power that needs 
to be continuously dissipated deep in the convective interior 
of hot Jupiters to explain their inflated radii (e.g., Gu, Boden- 
heimer & Lin 2004, Burrows et al. 2007, Ibgui et al. 2009). 
While Batygin & Stevenson (2010) emphasize such a deep de- 
position scenario, our own ohmic dissipation models say little 
about this scenario since we only calculate currents and ohmic 
dissipation in the superficial atmospheric region of the planet. 
However, another means by which hot Jupiter radii can be in- 
flated is by slowing down their rate of contraction, through 
modifications to the thermal structure of their overlaying at- 
mospheres which act as boundary conditions for the cooling 
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isentropic interiors. Indeed, Guillot & Showman (2002) ar- 
gue that, if ~ 1 % of the stellar insolation flux were deposited 
at pressures of tens of bars deep in the atmosphere, this would 
slow down cooling sufficiently to explain the inflated radii of 
hot Jupiters. While Guillot & Showman (2002) suggested that 
a downward flux of kinetic energy could in principle achieve 
such energy deposition, our Fig. 2 indicates that this is in fact 
naturally achieved by ohmic dissipation in our atmospheric 
models, with or without drag, for a magnetic field strength 
B>10G. 

More specifically, Guillot & Showman (2002) describe an 
evolutionary model in which 2.4 x 10 20 W are deposited at 
a location centered around 21 bar in the atmosphere of HD 
209458b, which can explain its observed inflated radius. This 
is achieved by the two ohmic dissipation models represented 
by the upper dashed and solid lines in our Fig. 2. In fact, since 
even more ohmic power is dissipated deeper in, we anticipate 
that models with fields even weaker than 10 G, possibly as low 
as 3 G, will be able to meet the inflated radius requirement for 
HD 209458b. This leads us to conclude that ohmic dissipa- 
tion deep in the atmospheres of hot Jupiters, which indirectly 
taps into the kinetic energy of dragged winds driven by stel- 
lar insolation higher up, is a promising scenario to explain the 
inflated radii of hot Jupiters. 

3. SUMMARY AND CONCLUSIONS 

In this Letter, we have computed the rate of ohmic dissipa- 
tion in representative, three-dimensional atmospheric circu- 
lation models of the hot Jupiter HD 209458b. We find that, 
for a fiducial magnetic field strength of 3 G, ohmic dissipa- 



tion starts dominating over stellar insolation heating at levels 
deeper than a few bars. The spatial non-uniformity of this ex- 
tra source of heat could induce Joule-driven circulation in the 
deep layers traditionally considered as "inert". For a magnetic 
field strength > 10 G, our models also indicate that enough 
heat is deposited at pressures of several tens of bars to slow 
down cooling sufficiently that the inflated radii of hot Jupiters 
can be explained. 

Our results hence suggest that magnetic interactions in hot 
Jupiter atmospheres play a fundamental coupling role for the 
dynamics and the thermal evolution of these planets. As such, 
our work calls for the problem to be treated more consistently: 
magnetic drag affects wind speeds and induces currents, while 
the ohmic dissipation of these currents can generate a deep 
atmospheric circulation which could, in turn, feedback on the 
circulation higher up. Indeed, the extra heat source might also 
enhance convection in the night side, which might in turn en- 
hance cooling. These various ingredients will have to be in- 
corporated consistently in circulation models for a better as- 
sessment of their influence on the structure and the evolution 
of hot Jupiters. Some diversity may naturally arise from vari- 
ations in the magnetic field strength and geometry of different 
planets. 



We thank Tristan Guillot for useful discussions, and Adam 
Burrows, Jeremy Goodman and the referee Douglas Lin for 
helpful comments on our manuscript. This work was sup- 
ported in part by the National Science Foundation under Grant 
No. PHY05-51164. 



REFERENCES 



Baraffe, I., Chabrier, G., Barman, T. 2010, Reports on Progress in Physics, 

Vol. 73, Issue 1, pp. 016901 
Batygin, K., Stevenson, D. J. 2010, ApJL, 714L, 238 

Burrows, A., Orton, G. to appear in "Exoplanets", Spring 2010 Space Science 
Series of the University of Arizona Press (Tucson, AZ); Ed. S. Seager; 
eprint arXiv: 0910.0248 

Burrows, A., Hubeny, I., Budaj, J., & Hubbard, W. B. 2003, ApJ, 594, 545 

Burrows, A., Rauscher, E., Spiegel, D. S., & Menou, K., eprint 

arXiv: 1005.0346 

Charbonneau, D., Brown, T., M., lathm, D. W., & Mayor, M. 2000, ApJ, 529, 
L45 

Christensen, U. R., Holzwarth, V. & Reiners, A. 2009, Nat., 457, 167 
Cooper, C. S. & Showman, A. P. 2005, ApJ, 629, L45 
Dobbs-Dixon, I., Lin, D. N. C. 2008, ApJ, 673, 513 
Draine, B. T., Roberge, W. G., Dalgarno, A. 1983, ApJ, 264, 485 



Hansen, B. M. S. 2008, ApJS, 179, 484 

Gu, P.-G., Bodenheimer, P. H. & Lin, D. N. C. 2004, ApJ 608, 1076 

Guillot, T., Showman, A. P. 2002, A&A, 385, 156 

Ibgui, L., Spiegel, D. S. & Burrows, A. 2009, eprint arXiv:0910.5928 

Iro, N, Bezard, B., & Guillot, T. 2005, A&A, 719, 727 

Laine, R. O., Lin, D. N. C, Dong, S. 2008, ApJ, 685, 521 

Liu, J., Goldreich, P. M., Stevenson, D. J. 2008, Icar., 653, 664 

Perna, R., Menou, K„ Rauscher, E., 2010, ApJ, 719, 1421 

Rauscher, E. Menou, K. 2010, ApJ, 714, 1334 

Sanchez-Lavega, A. 2004, ApJ, 609, L87 

Seager, S. & Deming, D. 2010, ARA&A in press, arXiv: 1005.4037 
Showman, A. P., Y-K. Cho, J., Menou, K. 2010, to appear in "Exoplanets", 

Spring 2010 Space Science Series of the University of Arizona Press 

(Tucson, AZ); Ed. S. Seager; eprint arXiv: 0911.3170. 
Showman, A. P. et al. 2009, ApJ 699, 564 



